Percolation-type description of the metal-insulator transition in two dimensions 
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A simple non-interacting-electron model, combining local quantum tunneling and global classical 
percolation (due to a finite dephasing time at low temperatures), is introduced to describe a metal- 
insulator transition in two dimensions. It is shown that many features of the experiments, such 
as the exponential dependence of the resistance on temperature on the metallic side, the linear 
dependence of the exponent on density, the e'^ /h scale of the critical resistance, the quenching of 
the metallic phase by a parallel magnetic field and the non-monotonic dependence of the critical 
density on a perpendicular magnetic field, can be naturally explained by the model. 
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The experimental observation of a metal-insulator 
transition in two dimensions has been a subject of 
extensive investigation, since it is in disagreement with 
the predictions of single-parameter scaling theory for 
noninteracting electrons Several theories, based on 
the treatment of disorder and electron-electron interac- 
tions by Finkelstein Q , have been put forward |^ . Other 
approaches considered spin-orbit scattering percola- 
tion of electron-hole liquid or scattering by impurities 
p| . To date there is no acceptable microscopic theory 
that describe quantitatively the observed data. 

Here we present a simple non-interacting electron 
model, combining local quantum tunneling and global 
classical percolation, to explain several features of the 
experimental observations. The main observations we 
want to understand are the following: 

• As the system is cooled down the resistance of sam- 
ples with density higher than some critical density 
extrapolates to a finite value at zero temperature, 
while that of samples with lower density diverges. 

• The resistance of the sample with the critical den- 
sity does not depend on temperature (at least for 
a limited range of low temperatures) . 

• The conductance of the critical-density sample is of 
the order of e^/h. 

• On the metallic side the functional dependence 
of the resistance is of the form R(T) = Rq + 
Ri exp{—A/T). The parameter A varies linearly 
with the density and vanishes at the transition. 

• In perpendicular magnetic fields this transition is 
continuously connected with the quantum Hall- 
Insulator transition jl^ . The critical density varies 
nonmonotonically with magnetic field, with a min- 
imum around v — 1. 



• Parallel magnetic fields destroy the metallic 
at least for densities near the transition [0,E 



)hase. 



Before introducing the model let us mention three other 
experimental observations - (a) The strong disorder is 
crucial to see the transition. In GaAs the transition is 
seen only at samples with low mobility (even with the 
same density). In fact, Ribeiro et al. have recently 



observed a zero-field metal-insulator transition in high- 
density n-type GaAs sample with strong enough disor- 
der (which was introduced by a matrix of randomly dis- 
tributed quantum dots), (b) There is additional exper- 
imental indications that the transition is not driven by 
interactions - Yaish and Sivan studied a system of 
two parallel gases, one of electrons and one of holes. The 
observed metal-insulator transition in the hole gas de- 
pended only slightly on the electron density, even though 
one expects that increasing electron density will screen 
the interactions between holes and suppress the metallic 
phase in the hole gas. On the contrary, increasing elec- 
tron density led to increasing conductance in the hole 
gas, indicating that its main role is to screen the impu- 
rity potentials in the hole gas. (c) There is a growing 
experimental evidence that even at the lowest available 
temperatures, the dephasing length, L^, is finite [p^|-p7|. 

Based on all these observations we now suggest the fol- 
lowing scenario to explain the experimental observations 
- the potential fluctuations due to the disorder define 
density puddles of size or larger in which the electron 
wavefunction totally dephases. (Density separation into 
puddles in gated GaAs was indeed observed experimen- 
tally by Eytan et al. using near-field spectroscopy.) 
Locally, between these puddles, transport is via quan- 
tum tunneling through saddle points, or quantum point- 
contacts (QPCs). Since between such tunneling events 
dephasing takes place, the conductance of the system will 
be determined by adding classically these quantum resis- 
tors. A related model was introduced by Shimshoni et 
al. to describe successfully transport in the quan- 
tum Hall (QH) regime. In fact, it was later concluded 
po| that the observation of a finite, quantized Hall re- 
sistance in the Hall insulator phase |^,^ can only oc- 
cur when the dephasing length is smaller than the size 
of these puddles - otherwise the Hall resistance diverges 
|p3| . In addition, this model also accounted for the ob- 
served current- voltage duality around the transition , 
a duality which was also observed in the zero-field tran- 
sition . The percolative nature of the system in the 
QH regime was indeed verified experimentally [|4j. Fi- 
nite dephasing length can also explain the observed non- 
critical behavior of the resistance near the QH-insulator 
transition ||2^. We will return to the QH regime below. 

We characterize each saddle point by its critical en- 
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ergy Cc, such that the transmission through it is given by 
T(e) = 0(e — Cc)- Thus the conductance through each 
QPC is given by the Landauer formula, 



h l + exp[(ec-/i)/fcT]' ^ ' 

where fj, is the chemical potential, and fpo is the Fermi- 
Dirac distribution function. 

The system is now composed of classical resistors, 
where the resistance of each one of them is given by (||) , 
with random QPC energies. In the numerical data pre- 
sented below, we solved a 20x20 system of QPCs (which, 
for simplicity, has the topology of a square lattice), each 
averaged over 1000 realizations of disorder, where the 
QPC energies were taken from a square distribution of 
width W. At zero temperature the conductors have ei- 
ther zero conductance or a conductance equal to 2e^//i 
and one has the usual second-order percolation transi- 
tion. The critical conductance exponent t is known at 
two dimensions and is equal to 1.3 M. In Fig. 1 we 
fit the experimental data of |Q] and of^to the expected 
critical dependence. Clearly, the agreement with the clas- 
sical percolation prediction is excellent. Moreover, the 
fact that while the density scale is so different between 
the two experiments, the conductance scale is identical, 
clearly demonstrates that e'^/h is the only conductance 
scale in the system. 
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Fig. 1. Comparison of the lowest temperature data of [g| (two 
sets of data, triangles and squares, 330mK, density given by 
the lower axis) and of Q (circles, 57mK, density given by 
the upper axis) to the prediction of percolation theory (solid 
line). Inset: Logarithmic derivative of the data which gives 
a line whose slope is the critical exponent. The percolation 
prediction (t = 1.3) is given by the solid line. For comparison 
a.t = \ slope is also shown (broken line). 

As temperature increases, the Fermi-Dirac distribution 
is broadened. Consequently the conductance of the insu- 
lating QPCs (ec > m) increases exponentially, while that 



of the transparent ones (ec < \i) decreases exponentially. 
Thus we expect to see rather dramatic effects as a func- 
tion of temperature. This is indeed depicted in Fig. 2. As 
temperature is lowered systems with slightly different re- 
sistance at high temperatures will diverge exponentially 
with decreasing temperatures. The resistance of systems 
on the metallic side (n > Uc) will saturate at zero tem- 
perature, while that of insulating samples will diverge. 
Note that there is an upward turn even on the metallic 
side of the transition. We will come back to this point 
below. The high-temperature resistance of the critical 
density network is naturally around h/e'^, the only resis- 
tance scale in this model. 
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Fig. 2. Temperature dependence of the resistance for systems 
of different densities. Below the critical line (bold curve) all 
curves saturate at zero temperature, while above it the resis- 
tance diverges. The resistance of the more metallic samples 
decreases exponentially (inset). 

For systems of exponentially distributed resistors the 
resistance of the whole circuit is determined by the crit- 
ical resistor - the worse resistor in the minimal perco- 
lating network |^^. Then the resistance of the network 
would be equal to the inverse of the Fermi-Dirac function, 
{1 + exp[(ec — fi)/kT]}h/e'^, where ec is its threshold en- 
ergy. Then clearly the overall resistance will be of the 
form observed experimentally, R = Rq + Ri exp(— A/T), 
with A varying linearly with the density and vanishing 
at the transition. In our case, the resistors on the metal- 
lic side have a bound distribution, and accordingly there 
will be other resistors, in parallel and in series, that will 
contribute to the overall resistance of the circuit. This 
will not change the functional form, but will renormalize 
the parameters i?o, and A. Such a functional depen- 
dence on the metallic side is indeed found numerically 
and displayed in the inset. In fact, close to the transi- 
tion, on the metallic side, as temperature increases, some 
QPCs that before had zero conductance, start to conduct 
and add to the overall conductance. Since the critical 



2 



percolation cluster is very ramified (in fact of fractal di- 
mension), there will be many such resistors in parallel to 
the main conducting network, and the effect of improv- 
ing these resistors will overcome the fact that resistors on 
the conducting network itself become worse. This leads 
to a downward turn of the resistance with increasing tem- 
perature even on the metallic side, the details of which 
may depend sensitively on the geometry. Only deeper 
into the metallic regime, as seen in Fig. 2, the overall 
resistance increases with increasing temperature. This 
also suggests that the density at which the resistance is 
approximately temperature independent is not the true 
critical point, but rather deeper on the metallic side. This 
is clearly seen in Fig. 3, where one can see a point where 
all the low-temperature curves nearly cross, well inside 
the metallic regime. The above discussion suggests that 
one should be cautious in associating the critical point 
with the "temperature- independent" point, as done rou- 
tinely in the experiment interpretations. 
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Fig. 3. Conductance vs. density (Fermi energy) for several 
temperatures. There is a density, well above the true critical 
point, where the curves seem to cross each other. 

We turn next to the effects of a magnetic field. The 
effect of a parallel field is straightforward to understand, 
as there have been several studies of transport through 
a QPC in parallel fields These experimental and 

theoretical studies demonstrated that the threshold den- 
sity where the QPC opens up increases parabolically with 
the in-plane magnetic field. This effect was attributed to 
coupling of the in-plane motion to the strong confinement 
in the vertical direction, leading to an increase in the en- 
ergy levels. Since this increase occurs for all QPCs in the 
sample, such a field in our case will strongly inhibit the 
metallic behavior: QPCs which were conducting at zero 
field, will have exponentially small conductance with in- 
creasing field. Once the density of conducting QPCs falls 
below the critical density, the system becomes an insula- 
tor, in agreement with experimental observations. 

The situation in perpendicular magnetic fields is more 
interesting, as QH states are formed. Transport through 



a single QPC in perpendicular field and the crossover be- 
tween the zero field limit and the QH limit have been 
studied in detail As expected, one finds that the 

critical energy oscillates with magnetic field due to the 
depopulation of Landau levels. In our case, we expect the 
oscillations to be smoothed out by the disorder and by 
the averaging over many QPCs. Thus only the strongest 
oscillation, near v = 1, may survive, leading to a single 
dip in the critical density vs. magnetic field plot, as was 
observed experimentally. This is indeed in agreement 
with our numerical calculation. We studied the energy 
levels of one puddle of electrons, which we modeled by a 
circular disk, in the presence of disorder In Fig. 4 
we plot the "critical density" - the number of electrons 
that need to occupy the puddle, so that the energy of the 
highest-energy electron will be enough to transverse the 
QPC equivalent in the bulk system to the critical 
density - as a function of magnetic field. Indeed we see a 
dip near = 1 with all other oscillations smoothed out by 
the disorder. This curve has a strong resemblance to the 
experimental data [T^ ] (inset). In addition, we expect 
that as the magnetic field is lowered below the ~ 1 
minimum, more than one channel will transverse some 
QPCs, leading to an increase in the critical conductance, 
as indeed reported experimentally. 
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Fig. 4. The critical density - the number of electrons in 
the puddle, so that the topmost energy will allow transport 
through the point contact - as a function of magnetic field, 
in the presence of a finite disorder. The continuous curve is 
an averaged fit through the (necessarily integer) data points. 
Inset: the corresponding experimental data 

All the above results and discussion demonstrated that 
many of the experimental observations can be explained 
in the context of the simple model introduced here. Nev- 
ertheless there are clearly other physical effects that need 
to be included in order to have a full picture of the ex- 
periments. In particular, electron-electron interactions 
are expected to play an important role in these low den- 
sities. As we can regard the metallic puddles described 
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above as quantum dots, one can use the abundant infor- 
mation about the role of interactions in such structures 
psf , to gain additional understanding of the character- 
istics of the puddles and the phase separation. Other 
effects, including the energy dependence of the trans- 
mission coefficient and the possibility of more than one 
channel through the QPCs, the temperature dependence 
of the dephasing length, the role of interband-scattering 
1^ and temperature-dependent impurities [|ll| may also 
be important to understand quantitative aspects of the 
data. Nevertheless, the fact that several important as- 
pects of the experimental data can be explained in the 
context of a simple model is quite encouraging. We ex- 
pect that the model presented here apply mostly to the 
GaAs samples. 

The picture described above can be checked exper- 
imentally. The experiments verifying the percolative 
structure in the QH regime |Q can be extended to the 
zero-field systems. (An experimental evidence for phase 
separation was observed at zero field in [Q.) An even 
more direct evidence of the percolative nature of the sys- 
tem will be local probes [Q. In fact, an enhancement 
in the fluctuations of the local chemical potential has al- 
ready been observed as the system enters the "insu- 
lating" phase, in a similar fashion to the enhancement of 
chemical potential fluctuations with closing of the barri- 
ers forming a quantum dot [Q. In fact, a "smoking gun" 
verification of the picture presented here, will be periodic 
oscillations of the local chemical potential on the insulat- 
ing side, due to depopulation of the Landau levels, as was 
observed in quantum dots |^ . 
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